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Abstract 

Recently an electric and magnetic field correlator appearing in the description of 
the heavy quarkonium system was evaluated on the lattice. Here, we give a non- 
perturbative analytical evaluation of this field correlator using a dual description of 
long distance Yang-Mills theory and using the stochastic vacuum model. The two 
predictions are both compatible with lattice data but show a different dependence 
on the quark separation. We discuss the analytic results in relation to the lattice 
data. 
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1 Introduction 



The semirelativistic interaction for a heavy quark-antiquark system is known 
in terms of electric and magnetic field correlators [1-7]. Due to confinement, 
these correlators cannot simply be evaluated in perturbation theory. However, 
it is possible to obtain them with lattice simulations. Analytic estimates of 
the correlator behaviour can be obtained by making some assumptions on the 
nature of the QCD vacuum or on the confinement mechanism. Until now the 
lattice results are compatible with, up to the available precision, the analytic 
calculations of the field correlators involved in the static and spin-dependent 
quark-antiquark semirelativistic interaction (see e.g. [4]). Recently, the correla- 
tors entering the velocity-dependent spin-independent part of the interaction 
were obtained on the lattice [8] and are compatible with the existing predic- 
tions (see [3,5-7]). 

In [6] it was shown that a spin independent term, - ( — ^ H = ] AV a , arises 

8 \mf wij J 

at order 1 / m 2 in the quark-antiquark potential from the Darwin and the spin- 
spin interaction. A previous analytic determination of AV a in the minimal 
area law model [6] gives zero. The first lattice QCD evaluation of AV a has just 
been carried out and gives non-zero negative value [8]. For quark-antiquark 
distances greater than 0.2 ~ 0.3 fm its behaviour is compatible with a negative 
constant with a value between and —1.5 GeV 3 . For distances shorter than 
0.2 fm AV a falls off like 1/r. 

In this letter we calculate AV a using a dual description of long distance Yang- 
Mills theory [9]. We give also the predictions of the stochastic vacuum model 
[3] thereby completing a recently given comparison between these models [4]. 
Finally we compare the theoretical predictions with the lattice measurements 
and discuss the implication of these results for the confinement mechanism. 



AV a is defined as following 



1 . 



AV a = AV a B + AV a E , (1) 

2 T/2 T/2 



A V s = — lim % 



J dt J dt'((B(z 1 ) ■ ~B(z[))) - ((B(zi))) • ((B(z[))) , 

-T/2 -T/2 

(2) 



Here and in the following e stands for the usual QCD coupling constant, q s = 
e 2 /4-7r, while g denotes the dual coupling constant g = 2w/e. We have used this 
unconventional notation to make contact with the dual theory papers. 
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T/2 T/2 



AVf = lim 



(ft / ^((E(^)-E(4)))-((E(^)))-((E(^))), 



-T/2 -T/2 



where 



(3) 



^ = F 0i , fii = F 0i , = <9^A, - d v Afj, + ie [A^ A v \ , 

piiv _ \ Axvpa rp 

1 2 pCT ' 



and 



«/(^)» 



J VAe iSYM(A) Ti P < 


/(A) exp 


ie <j) dx^ 1 A^(x) 
. r 


1 


y PAe i5YM(A) TrP exp 


L r 


f dx^A^x) 





(4) 



We define, also, the Wilson loop: 



W(T) = P exp 



ie j> dx^A^x) 
r 



(5) 



The closed loop T is defined by the quark (anti-quark) trajectories z 1 (t) (z 2 (£)) 
where t varies from the initial time —T/2 to the final time T/2. The quark 
(anti-quark) trajectories Zi(t) (z 2 (t)) define the world lines F\ (r 2 )- The world 
lines Ti and T 2 , along with two straight-lines at fixed time connecting the end 
points, make up the contour T. We use also the notation z\ = z\{t) and 
z[ = zi(t'). As usual A M (x) = A'j I (x)\ a /2, Tr means the trace over colour 
indices, P prescribes the ordering of the colour matrices according to the 
direction fixed on the loop and Sym(A) is the Yang Mills action including a 
gauge fixing term. 



AV a is a gauge invariant quantity. 



The field-strength expectation values in (2), (3) can be expressed as functional 
derivatives of log(W(r)) with respect to the path T [5]. In fact the change 
induced in (W(r)) by letting z{{t) -> z£(t) + 6z£(t) where 6z?(-T/2) = 
5z?(T/2) = is 

5i\og(W(T)) 

(fc(*i)» = eSSia , {zi) » ( 6 ) 



2 



5S' IU (z 1 ) = {dz^8z\ - dz^Sz?) , (7) 
and, varying the first quark line again, we have 



Hence, the evaluation of (2) and (3) follows straightforwardly from the as- 
sumptions on the Wilson loop behaviour. 

In the definition of AV a (eqs (2), (3)) the quarks are static, i. e. z x and z 2 are 
fixed in time. Accordingly the lattice calculations have been done inserting 
field-strength tensors in a rectangular Wilson loop. 

In the static limit we note the following identities involving the static potential 
[6]: 



T/2 T/2 



J dtVV = e J dt ((E(zi))) , (9) 

-T/2 -T/2 
T/2 T/2 T/2 

J dtAV = -te 2 J dt J dt'((E(z 1 )-E(z' 1 )))-((E(z 1 )))-((E(z[))) 

-T/2 -T/2 -T/2 

T/2 

+e J dt((D-E( Zl ))), (10) 

-T/2 

where _D M is the covariant derivative. Therefore, we have the following alter- 
native way to express AVj 3 and AVj^: 



T/2 T/2 



/ * / w 



-T/2 -T/2 
T/2 T/2 

, tp , i r '( 5 2 log(H/(r)) 
AVF = - lim - dt dt' 



1 " t->co T J J SSoiizJSSoiW 

-T/2 -T/2 
T/2 T/2 



- J dt AVo + e J dt ({D ■ E{ Zl ))) . (12) 



-T/2 -T/2 

The last term in (12) is the Darwin potential [5,7]. 
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2 Electric-magnetic duality 



We assume that for large loops the Wilson loop W is the same as the Wilson 
loop W e fi determined by an effective dual theory which is weakly coupled at 
long distances [7]: 

(W(T)} = (W cS (T)} (13) 



where 



J VC VB exp 


i J dxC eS {G^ v ) + £ G f 




J VC VB exp 


i 


[ dxjC eS (G% = 0) + £ GF 



(14) 



£ e fr is the effective Lagrangian describing the interactions of dual potentials 
C M and monopole fields B^ the dual potentials interact with quarks via the 
Dirac string tensor G^ u : 

/f dy^ dy 
ds J dr -q^-q^ S ( x ~ y( s ' r )) > ( 15 ) 



where y(s, r) is a world sheet with boundary V swept out by the Dirac string. 
The colour electric field D and the color magnetic field H are the compo- 
nents of the dual field-strength tensor Q^ u = d^C u — d u C^ + G^ : Di = Cj i 
and Hi = Q 0i . The Duality assumption (13) yields expressions for the Yang 
Mills expectation values (6) and (8) in terms of expectation values of the 
corresponding operators in the dual theory [7] 



'<B(^» = ^((H^eff , (16) 

((E^H^D^eff, (17) 
((B( Zl ) ■ B(zJ)» - ((B( Zl ))) • «B(4)» = -\i ^is^ , (18) 

((E( Zl ) ■ E(^)» - «E(^)» • «E(^)» = - |i H f^ eS , (19) 

where Df = G^, Hf = G^. In the classical approximation to the dual theory 
(justified because the theory is weakly coupled at large distances) all quantities 
on the right-hand side of equations (14), (16)-(19) can be expressed in terms 
of the solutions of the classical equations of motion, in particular 
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i log(W(r))eff = - 1 dxC ef[ (G%) , (20) 

«D(x)» eff = D(x), (21) 
«H(*)» cfr = H(x). (22) 

The explicit evaluation of (18) in the classical configuration gives [7]: 

AVf = -\e 2 V x • V x ,G NP (x,x')| x=x , =zi . (23) 

where the function G NP satisfies the equation 

(-A x + 6^ 2 (x)) G NP (x,x') = - ff^, , (24) 

v ' 47T | X — X 

with g = 2ir/e and B is the static monopole field. 

Now we use equation (24) to obtain the behaviour of G NP near the quark 
sources. Let us introduce the vector 



K(x) = 



NP/ „/\ 



(25) 



In terms of cylindrical coordinate with origin in z±, we have x— zi = z e z +pe p , 
where the unit vector e z is directed along r = zi — Z2 and the unit vector e p 
is orthogonal to e z in the plane of x and r. By symmetry 

K = K z e z + K p e p . 

As a consequence, (24) can be splitted into two equations, one for each com- 
ponent: 

(-A + 6^( X ))^-g^g, (26) 

(-A + «^<x))*, = -*g^, (27, 

1 d d 2 1 d 2 d 2 ~ 1 
where, now, A = -— + — - + -^tttt + ttt and A = A -. Let 6 be the 

p Op Op z p z 0(f) z oz A p A 

angle between r and x — Zi, where for # = 7r the vector x — zi is along r, so 
that z = |x — zi| cos#, p = |x — z± \ sin 9. On the string, B must vanish, and 
its behaviour near quark 1 is given by [10] 

i3 = 6v / l + cos^|x-z 1 | (v ^- 1)/2 , (28) 
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where b is a function weakly dependent on r (the quark-antiquark distance). 
In Table 1 we list values of b obtained by numerical integration of the static 
equations for B for varying separations r. Putting (28) in (26)- (27) and keeping 
only terms which are singular near the quark 1 (for x — > Zi), we obtain 



_ A^ = -^Acos^(l + cos^)|x- Zl |^- 3 , 
2n 

- AK P = - sin 6 ( 1 + cos 9) |x - Zl | . 
2n 



(29) 
(30) 



The solutions of these equations are [10] 



3g 2 b 2 
2tc 
3g 2 b 2 
2vr 



K z 



cos 9 cos 2 9 
+ 



3 l-y/3 3 + V3 
1 cos 9 



^ sin^lx-zil^ 1 . 



(31) 
(32) 



Finally 



AVf = -^ 2 V x .V x ,G NP (x,x') 



X=X'=Zl 



,-e 2 V x • K(x) 



X=Zl 



= -8 7r6 2 (v / 3 + 2)|x-z 1 |^- 2 



(33) 



X=Zl 



The result (33) is divergent. However the dual theory has a natural ultraviolet 
cut-off which is the mass M rs 0.640 GeV of the dual gluon, and we can put 



x — Zi 



V3-2 



M 



(34) 



X=Zl 



From this prescription and the data of Table 1, we can derive the values of 
AV a B at different distances as given in Table 2. We observe that in the range 
of interest (0.2 fm < r < 1 fm) the magnetic contribution to AV a can be 
considered almost a constant —0.3 GeV 3 ). 

In order to evaluate AV^ in the dual theory we use the following analogue of 
eq. (10) in the dual theory: 



lim 

T 



1 

oo T 



T/2 T/2 



dt 



-T/2 -T/2 



4 MDfa))) 
3 SDHz[) 



cS 
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T/2 T/2 

/ rftA\/ + ^e | dtV ■ ((-D( Zl ))} cS . 

-T/2 -T/2 



(35) 



Comparing (35) with (12) we obtain 



T/2 T/2 



J dt ((D • E(zi))) = - / dtV-((D(^)» 



eff • 



(36) 




Eq. (36) expresses the Darwin potential in terms of quantities calculated in 
the dual theory and has a simple interpretation. The left hand side is the 
expectation value of the electric current operator in Yang Mills theory and 
the right hand side is the corresponding quantity in the dual theory (note 
that (36) and (17) are independent results valid only on the quark trajectory). 
Since in the dual theory the colour electric field D is given byD = VxC + D 5 
the dual potential does not contribute to divergence of D and we have 



where r = z± — Z2. The dynamics of the dual theory fixes VxD, that is, the 
monopole current which, because of (36), does not contribute to the Darwin 
term. Thus the monopole degrees of freedom, which in the classical approxi- 
mation to the dual theory produce the dual Meissner effect and confinement, 
contribute to AV A only via the first term on the right hand side of eq.(35). 
Thus comparison of AV a predictions with lattice data provides a crucial test 
of the dual description of long distance Yang Mills theory. 

Putting (37) in (35) we obtain 



where a s = e 2 /Air. The quantity G NP determining AVj 3 (33) is just the nonper- 
turbative part of the monopole potential and so (33) is the natural magnetic 
counterpart of (38). Thus both the electric and magnetic parts of AV a directly 
measure the characteristic feature of the dual picture. The static potential V 
obtained by solving numerically the equations of motion can be represented 
by the following analytic parametrization [11]: 



V-((D(^))) cff = V-D(^) = e5(r), 



(37) 



AVf = -AV + ^a s <y(r) = -A( 




) 



(38) 



V (r) 



4«s 
3 r 




+ or — 0.646 v/cra! 



(39) 
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where a rs 0.2 GeV 2 is the string tension. Substituting (39) in (38) we obtain 
AV* = 1 (0.511) 2 V 0511 -2°-. (40) 

In Table 3 we list some values of AVf for different quark-antiquark distances. 
Notice that the dependence on r is strong. In particular the short-range be- 
haviour (r — > 0) is like \jr. In Table 4 we sum up the magnetic and electric 
contributions and give the complete AF a behaviour. In the long-range region 
(r — > oo) A14 has only magnetic contributions with an asymptotic value of 
0.2 ~ 0.3 GeV 3 . In the short-range region AV a has a behaviour like 1/r. 



3 Stochastic vacuum model 

In the stochastic vacuum model, one assumes that, in Euclidean space, the 
Wilson loop has the following behaviour [3] : 

\og(W(T)) = 

-\ J dS pu (u) J dS Xp (v)(5, x 5 up - 5 pp 5 vX ) (D(h 2 ) + D^h 2 )) 

s s 

d 

+ (h p h\5 vp - h p h p 5 uX + h u h p S pX - Kh x 5 pp ) —D^h 2 ) , (41) 

where h = u — v. The functions D and Di, called correlator functions, can be 
interpreted as vacuum form factors. The perturbative contributions (contained 
in Di) can be evaluated by the usual e 2 expansion of the Wilson loop, while 
the nonperturbative contributions, dominant in the long-range behaviour come 
only from lattice calculations [12]. The integration is performed over the mini- 
mal area surface S with contour T. Up to order 1/m 2 the minimal surface can 
be identified exactly with the surface spanned by the straight-line joining the 
point (t,zi(t)) on the quark line 1 with the point (t, %2(t)) on the quark line 
2 of T. Therefore, for our purposes, the generic point w M (t, s) on the surface S 
will be given by: 

u° = t u = szi(i) + (l-s)z 2 (f) (42) 

and the surface element dS pu {u) by 

dSij{u) = dt ds rj(t) , (43) 
dS lj (u) = dtds (sz u (t) + (1 - s) z 2i (t)) rj (t) , (44) 
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where t G [-T/2,T/2], s G [0, 1] and r = z x - z 2 . 
From equation (8) and (41), expanding 

Zl (t') = z l (t) + (t-t')z l (t) + --- , (45) 

we have in the static limit (i. e. neglecting the velocity dependent terms) 

#TWTM = "(Mm. - SuSmr) (D(1?) + Z>i(* 2 )) . (46) 

dbi m (zi)db rs (z' 1 ) v ' 

Putting this expression in (11) we obtain 

oo 

AV a B = -e jdt (D{t 2 ) + D 1 (t 2 )) . (47) 
o 

In an analogous way, since in the static limit 

sio g{ w(v) } _ _ 3 ( D{e) + Di(f) + e A. Dm ) t m 



SSoifa) SSoiiz'i) \ - ■ dt 2 

we obtain, from the first equality in equation (12) 

oo 

AV? = -6 f dt (Dp) + -Dtfj) , (49) 
o 

where a partial integration was performed. 

In the usual long-range parameterization, suggested by the direct lattice eval- 
uation of the correlator functions D and D x [4,12], 

D{x 2 )=de- Slxl , 5^1 GeV, d « 0.073 GeV 4 , (50) 

D 1 (x 2 )=d 1 e- Sllxl , <JiwlGeV, d x 0.0254 GeV 4 , (51) 

the electric and magnetic contribution to A14 are 

AK B = -6^ - 6^ w -0.59 GeV 3 , (52) 

AVf = -6^ - 3$- « -0.51 GeV 3 , (53) 
5 5i v ' 
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and then 



AK = AV® + AV a E « -1.1 GeV 3 



(54) 



In contrast with the dual theory the long range contributions to AV a in the 
stochastic vacuum model are self energy corrections independent of r. For this 
reason these terms were neglected in [4] and did not appear in the minimal 
area law evaluation of [6]. We expect that a r dependence will arise in the 
short range region, taking into account higher order perturbative contribu- 
tions. Since the r dependence is very weak also in the dual theory for the 
magnetic contribution AV a B , the difference between the two models is in the 
electric contributions AV a E and is mainly due to the Darwin term. In fact, in 
the static limit the dual theory predicts 

T/2 T/2 

e J dt ((D • E(zi))) = ^e 2 f dt5{r). (55) 

-T/2 -T/2 



Eq.(55) reflects the fact that in the dual theory the color electric field is 
generated by the monopole current which does not contribute to the right hand 
side of eq.(36). This contrasts with the corresponding result of the stochastic 
vacuum model: 

T/2 T/2 

e J dt ((D • E(zi))) = J dt VVo + self-energy contributions. (56) 

-T/2 -T/2 



4 Conclusion 



The dual theory predictions of AV a are in very good agreement with the lattice 
data. The asymptotic value for large distances is negative and very well inside 
the error bars of the lattice prediction. Surprisingly also the short range predic- 
tion, which in principle could be beyond the validity range of the dual theory, 
reproduces the 1/r behaviour of the lattice fit. Also the long range result of 
the stochastic vacuum model is negative and compatible with the lattice data, 
but of self-energy type. A r dependence arises from the perturbative contribu- 
tions only at the next-to-leading order. However the main contribution to the 
- behavior of V E in the dual theory comes from the second term on the right 
hand side of eq. (40), which is clearly nonperturbative. Furthermore, as is seen 
from Table 2, V a B also increases with decreasing r. After this calculation the 
main difference in the two models seems to be clear. In the dual theory only 
the quark charges contribute in the static limit to the Darwin potential, while 
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in the stochastic vacuum model the Darwin potential is related to the full 
static potential. This means that a truly accurate lattice determination of the 
field correlator (which is possible for the magnetic part of the correlator and 
it is already in progress [13]) is of great interest and could provide significant 
information concerning the nonperturbative quark interaction. 
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Table 1 

b vs r, the qq distance. The table shows a slight dependence. 



r (fm) 


b 2/ M Vs-i (QeV 2 ) 


0.175 


0.007 


0.26 


0.006 


0.35 


0.0055 


0.7 


0.0048 


1.4 


0.0045 



Table 2 

AVj 3 at the ultraviolet cut-off M = 0.640 Gev vs r, the qq distance (a s = 0.37 and 
cr = 0.2 Gev 2 ). 



r (fm) 


AV a B (GeV 3 ) 


0.175 


- 0.407 


0.26 


- 0.350 


0.35 


- 0.320 


0.7 


- 0.277 


1.4 


- 0.263 
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Table 3 

AV^ vs r, the qq distance (a B = 0.37 and a = 0.2 Gev 2 ). 



r (fm) 


AV* (GeV 3 ) 


0.175 


- 0.402 


0.26 


- 0.275 


0.35 


- 0.208 


0.7 


- 0.109 


1.4 


- 0.057 



Table 4 

AV a = AV^ + AVf vs r, the qq distance (a s = 0.37 and a = 0.2 Gev 2 ). 



r (fm) 


AV a (GeV 3 ) 


0.175 


- 0.809 


0.26 


- 0.625 


0.35 


- 0.528 


0.7 


- 0.386 


1.4 


- 0.320 
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